String Picture of a Frustrated Quantum Magnet and Dimer Model 
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We map a geometrically frustrated Ising system with transversal field generated quantum dynam- 
ics to a strongly anisotropic lattice of non-crossing elastic strings. The combined effect of frustration, 
quantum and thermal spin fluctuations is explained in terms of a competition between intrinsic lat- 
tice pinning of strings and topological defects in the lattice. From this picture we obtain analytic 
results for correlations and the phase diagram which agree nicely with recent simulations. 
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The interplay of quantum and thermal fluctuations in 
geometrically frustrated magnets raise new and fascinat- 
ing issues both from a theoretical and experimental per- 
spective 0, ■ Frustration prohibits the formation of a 
collinear Neel state, inducing a macroscopic degeneracy 
of the classical ground-state. Even for magnets with a 
discrete Ising symmetry, the complexity of the ground- 
state manifold may endow the system with a continuous 
symmetry. Such symmetry is of particular importance 
to 2D quantum magnets at low temperatures since then 
the Mermin- Wagner theorem applies, precluding an or- 
dered phase 3]. Another possible but contrary scenario 
is "order-from-disorder" where zero-point fluctuations se- 
lect a small particularly susceptible class of the ground- 
state manifold and yield an ordered symmetry-broken 
state 0. This poses the question if competing fluctu- 
ations about the classical ground-states can lead to new 
strongly correlated states and (quantum) phase transi- 
tions of unexpected universality classes. 

The possibly simplest realization of classical frustra- 
tion is provided by the antiferromagnetic Ising model 
on a triangular lattice (TIAF). This model is disordered 
even at zero temperature with a finite entropy density 
and algebraic decaying spin correlations ||. Quantum 
dynamics arise if a magnetic field is applied transverse 
to the spin coupling. For this transverse field TIAF, 
and its companions on other 2D lattices, Moessner and 
Sondhi have argued the existence of both ordered and 
spin liquid phases 0, Q- Experimental realizations of 
these models can be found either directly in magnets with 
strong anisotropy, e.g., in L1H0F4 Q, or indirectly (via 
a related (2+l)D classical model) in stacked triangular 
lattice antiferromagnets [U with strong couplings along 
the stack as studied in recent experiments on CsCoBr3 
. Transverse field Ising models also describe the singlet 
sector below the spin gap of frustrated antiferromagnetic 
quantum Heisenberg models [1(1 ] and phases of quantum 
dimer models (QDM) on the dual lattice whose order- 
ing depends on the lattice structure [Tlj . Interestingly, 
a disordered QDM appears to be a promising candidate 
for quantum computing fl2j |. 

In this Letter we explain the effect of simple quantum 
dynamics on a frustrated Ising model by mapping the 



transversal field TIAF to a stack of 2D lattices of non- 
crossing elastic strings. We show that the strings are de- 
scribed by a frustrated 3D XY model with a 6-fold clock 
term. Since our mapping yields the XY coupling con- 
stants, we can compute explicitly the phase diagram at 
arbitrary transverse field strength. The diagram is shown 
in Fig. Previous studies of the model started from 
a Landau- Ginzburg- Wilson (LGW) theory 0], which 
neglects frustration 14!], to describe the vicinity of the 
quantum critical point at strong transverse fields U If 
the phase of the complex LGW order parameter is iden- 
tified with the string displacement, the LGW and string 
actions agree except for frustration. However, the string 
approach shows that the configuration space is restricted 
since strings cannot cross each other which constrains 
their displacement. Moreover, the LGW approach does 
not provide the coupling constants, including the sign 
of the clock-term which is important for the selection of 
the ordered state. Very recent Monte Carlo simulations 
support the LGW prediction for the phase diagram but 
the computations were performed for the related (2+l)D 
classical model whose equivalence to the quantum prob- 
lem involves a singular scaling limit which makes simula- 
tions difficult 0|. Recently, the weak field behavior has 
been studied in terms of a quantum kink crystal jlf^ . 

Model for magnet, dimers, and strings. — The trans- 
verse field TIAF has the Hamiltonian 



(i,j) 



(i) 



with nearest neighbor coupling J, transverse field T, and 
Pauli operators a x , a z . First, we will consider the clas- 
sical ground states for T = which have one frustrated 
bond per triangle. A complete dimer covering of the dual 
hexagonal lattice is obtained if a dimer is placed across 
each frustrated bond 

G3, see Fig.EDt a). For each dimer 
configuration a height profile h can be defined on the 
triangular lattice sites as follows. Starting at an arbi- 
trary site with some integer number, one follows a trian- 
gle pointing down clockwise and changes the height by 
+2 (-1) if a (no) dimer is crossed. Repeating this process 
for all down pointing triangles, one obtains a consistent 
height on all sites, see Fig. ^d). The string representa- 
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FIG. 1: (color online) (a) Relation between spins and dimers 
(b) Vortex-anti-vortex pair, (c) The two flat states with flip- 
pable plaquettes (grey), (d) Mapping of dimers to strings 
where the numbers denote the height profiles. The string dis- 
placement (arrows) is determined by the height of the two 
plaquettes (grey) which are joined by the displaced dimer. 

tion follows from a given dimer state by the subtraction 
of a fixed reference state with all vertical bonds occupied 
[13, see Fig. Did). The strings fluctuate but remain di- 
rected (along the reference direction) and non-crossing 
(due to the frustration of the spin system) . They act as 
domain walls across which the height (after subtraction 
of the reference height) changes by 3. A shift of h by 3 
corresponds to a string translation by a mean string sep- 
aration t — 3a/2 with a the triangular lattice constant. 
Thus the string displacement is u — £h/3 + UQ with h the 
coincident (original) height of two hexagonal plaquettes 
which are joined by the displaced (non-vertical) string 
segment, cf. Fig. (Hd). uq is a global constant which 
differentiates between two non-equivalent states of flat 
strings that are not related by shifts by a lattice vector 
of the triangular lattice, see Fig. dc). For uq — na/2, 
n £ Z, the spins have orientations (+ H — ) on the three 
sublattices and the non-vertical dimers sit on straight 
lines that are locked between the triangular lattice sites. 

uq = a/4 + na/2 selects orientations (H 0) with the 

straight lines locked on the sites, i.e., one sublattice has 
alternating spins as indicated in Fig. ^c) by plaquettes 
with a zero. 

By coarse-graining one obtains a continuous field u(r) 
which allows to write the effective free energy of long- 
wavelength fluctuations of the string lattice in the form 
of a continuum elastic energy, 

T A = J d 2 r{^-(d x u) 2 + C f(d y uf + V L (u)} (2) 

with compression (en) and tilt (C44) modulus which are 
both of entropic origin. Vl(u) is a periodic potential 
which reflects the discreteness of the lattice. Since equiv- 
alent flat states are related by shifts of all straight lines 
by a/2 and u = ah/2 + uq, the locking potential must 



favor integer h, 

Vl = — vcos(2nh) 



-vcos (4tt(u — uo)/a) (3) 



with v > [18|. For Jcyyc^, > 2-n/a 2 the potential Vl 
is relevant (under renormalization) and the strings lock 
into that state which maximizes the entropy by allowing 
for local fluctuations that do not push the strings out of 
a tube of width I around their straight reference posi- 
tion, i.e., |u| < 1/2. This condition restricts spin flips 
to certain plaquettes [grey in Fig. [Jc)]. Taking into ac- 
count that strings are directed and non-intersecting, we 
find for a lattice of N sites that the states of orientation 
(+-0) and (+ + -) allow for 2 N / 3+N / 12 = exp(0.2888 N) 
and 18 N / 12 x 2( 13 / 18 ) 3Ar / 6 = exp(0.2844 N) localized con- 
figurations, respectively. Thus (H — 0) is selected and 
m = -a/4 in Eq.©. 

The classical spin correlations can be easily obtained 
within the string picture. They follow from the local 
string density p(r) = p(l — d x u(r)) with p — 1/t as 
{pi<Jj) — {e* p( T ^ dx y The compression modulus in 
Eq. (0) can be obtained from the equivalence of non- 
crossing strings to ID free Fermions QJI which is based 
on the Pauli principle and yields en = ir 2 p 3 /g with 
g = pen the string tension. Then the periodic poten- 
tial Vl is irrelevant since ^Jc~\\c^l — np 2 < 2ir /a 2 . g can 
be obtained from a random walk argument, however, its 
actual value is unimportant since the mean squared dis- 
placement {[u(r) — u(0)] 2 ) = 1/(7^011044) lnr becomes 
independent of g. This yields the known exact result |20| 
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String model. — Now we consider a finite transversal 
field. First, we use the Suzuki-Trotter theorem to express 
the partition function of Eq. in terms of a (2 + 1)D 
classical Ising model with reduced Hamiltonian [21| 
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K±aikaik+i (5) 



which consists of n antiferromagnetic triangular lattices 
with if || = J/(nT) which are coupled ferromagnetically 
with K± = |ln(nT/r). The correspondence becomes 
exact for an infinite Trotter number n — > 00. Next, we 
relate the classical spins aik = ±1 to the height pro- 
file hik which is defined for each layer as above. The 
starting height is fixed at one site in each layer by set- 
ting tiQk = 0, (3) if (Tofe = +1, (—1) along a column 
with layer index k. Then the height on all sites is fixed 
(modulo 6) for a given spin state since parallel (anti- 
parallel) spins imply a change of h by +2 (—1) from site 
to site if the down pointing triangles are traversed clock- 
wise. If one spin is flipped, h changes on that site by ±3. 
One easily proofs that these rules are met by the relation 
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a ik = cos(QR ife + 7r^ fc /3) = ±1 with Q = (47r/(3a),0,0) 
and lattice sites Rjfc. The intra-layer coupling maps to 
— cos[^(/ijfc — hjk with a shift rj^j = +1, —1 for in- 

plane bond directions (a, 0), (a/2, ±V3a/2), respectively. 
For the inter-layer coupling rjij =0. In each layer the 
spin states map again to a string lattice which, however, 
now contains vortex-anti-vortex pairs formed by trian- 
gles with 3 frustrated bonds which span string loops, see 
Fig. [IJb). Notice that the strings remain non-crossing 
since a triangle can have either 1 or 3 frustrated bonds. 
Including the lock- in potential of Eq.© with u — —a/4 



and since h = 3(u 
string Hamiltonian 



uq)/£, we obtain the reduced 3D 



Ho = -K 



<ij>,k 
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(u lk - Ujk + rjija/2) 



Uik — Uik+l, 



^y^cos 

i.k 



6tt 



(6) 



with v > and the shift rjij reflecting frustration. Since 
Uik can vary over the bonds of a triangle only by +a or 
—a/2, the energy is minimized for a non-uniform change 
of along the triangles which is distinct from a con- 
tinuous field whose oriented uniform change defines a 
helicity, giving rise to a Z2 symmetry j2^|. Thus the 
Ising model of Eq.© maps to a stack of 2D lattices of 
non-crossing strings which is described by a (2+l)D frus- 
trated XY model with a 6-fold clock term. This resem- 
bles the GLW theory if the string displacement is 
identified with the phase <p of the order parameter via 
(j) = 7tu/£. However, there are two important differences, 
(i) the in-plane XY coupling is frustrated and (ii) there is 
a topological constraint on <f> since u is restricted by the 
non-crossing condition which increases the phase stiffness 
on large length scales. The XY coupling allows for vortex 
loops which are in general superpositions of two types. 
Vortex-anti- vortex pairs [cf. Fig.^b)] occurring in many 
layers form loops oriented perpendicular to the planes, 
while the boundaries of 2D regions along which strings 
in adjacent planes are shifted by 21 form parallel loops. 
If the loop size is bound, the XY couplings of Eq. © 
can be expanded in Uik, and each layer is described by 
Eq. © with a harmonic interlayer coupling which can 
render the lock-in potential relevant. 

Phase diagram. — First, we note that the layers can- 
not decouple independently from the vortex unbinding 
transition in the layers since parallel loops of diverg- 
ing size can occur only at a critical value for the ef- 
fective if || which is by factor of 1/8 below that for the 
in-plane dissociation of vortices 23] . In addition, a p- 
fold clock term is irrelevant under renormalization if 
p > 3.4 |24|. and thus we expect a quantum phase tran- 
sition in the 3D XY universality class. To locate the 
transition, we factorize the partition function of Eq. © 
into a spin wave and a vortex part by the substitution 
jT[cos(0)-l] -» Em=-oo ex P[-^(^- 27rm ) 2 /2], denoted 



as Villain coupling 25] . The couplings K are known in 
two limits, K = K for K -> 00 and K = 1/(2 ki(2/K)) 
for K — * 0. In fact, these limits are realized for K± and 
Ku at large Trotter numbers n, leading to 



K 



iln(nT/r), 



K 



l -\n-\2nT/J). 



(7) 



At T = 0, one can set T ~ J/n so that for n — > 00 the 
coupling K± remains finite, and the system shows 3D 
behavior. For any finite T, however, K±_ must diverge 
with n and the system behaves effectively as 2D on large 
length scales. Eliminating n from Eq. (J7J, we obtain a 
relation between K\\ and K±_ which makes the parameter 
space ID. At large n, we can expand this relation, giving 



K\\K±= 2 [l-ln(2r/J)tf||] . 



(8) 



For layered XY models, dimensional crossover scaling |26j 
can be used to obtain a relation between the 3D critical 
value for the in-plane coupling and the corresponding 
for a system of n layers, 



1 K c 

n KS, 
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K \\ 



1/2 



1 - 



K c 
K c 



(9) 



with the 3D XY exponent v w 2/3 and a constant 7. 
This yields for n — 1 in the limit <C the following 
relation between and Kf, 



= K 



(10) 



At the quantum critical point the classical string system 
of Eq.© is 3D and thus at Ku = the relations of 
Eqs. (|HJ) and I|1U|) must be identical for consistency. This 
implies 7 = l/(2Kf) and T c /J = \e v l K ^. The usual 
Kosterlitz-Thouless (KT) argument yields for the vortex 
unbinding transition K\ = 2/n x 2/V3 on the triangular 
lattice [22J, leading to T c /J — 1.24. However, this es- 
timate neglects renormalization effects due to the clock 
term, frustration and the non-crossing of strings which 
should provide a net increase of r c . Interestingly, recent 
Monte Carlo studies suggest T c / J w 1.65 ||. For T < T c 
the clock term is relevant, and the system is ordered. 

At finite temperature, 2D XY physics dominate at 
large scales, and KT singularities are expected [2f|. Us- 
ing from Eq. (jTU|> in Eq.© with Ku — K° we obtain 
the critical coupling K^(K j_) as function of K± and large 
but finite n. This result can be used to determine the vor- 
tex unbinding transition at finite T. We set Co = e 2K± 
so that n = CoT/T in K^(K±). The renormalization 
of the coupling constants K±, Ku must be dependent, 
and from Eq. © the effective Kf = 1/(2 m(2Cl7J)) 

where £ = e 2K± is the renormalized Co ■ Now we can use 



Kl = Kf and K 



Kf and expand K^(K±) for large 
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T 2 from scaling i — e- 
J 2 from ti h — •' 
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FIG. 2: (color online) Phase diagram predicted by Eq. 1111 
and Monte Carlo results of Fig. 1 in Ref. A critical phase 
(C) is separated by a boundary at T c ,i from the paramagnetic 
phase (PM) and at T c ,2 = 4/9 T Ci \ from an ordered phase (O). 

£ which yields in the limit n — > oo the phase boundary 
of a critical phase (C) with bound defects, 

where b = je v / Kl £o (2 In C) -5 ^ 3 is a numerical constant 
which is fixed by the (unknown) renormalization of K± 
and remains finite for n — * oo. However, for consis- 
tency, since Co — nT/T the effective coupling must be- 
have at large n like K*± ~ (n/n c ) 3 / 5 with a charac- 
teristic n c = bK£(J/T)(T/r c ) which is a measure for 
the strength of quantum fluctuations and characterizes 
the effective system size along the Trotter ("imaginary 
time") axis. At a lower T Cj 2 there is a second transition 
to an ordered state where the clock-term in Eq. © be- 
comes relevant and locks the strings to the lattice. For 
a single layer, this transition is again of KT type with 
critical Kf = 9/(2tt) x 2/\/3 HjJ, and the boundary of 
the ordered phase can be obtained analogously to T C]1 
but with 7 replaced by jKf/Kf in b below Eq. (|llfl . 
yielding T Cj2 = (4/9)T Cj i. Close to the quantum critical 
point both temperatures vanish ~ (r c — r)" as expected 
from scaling. Fig. [21 compares the result of Eq. to 
recent Monte Carlo data for the phase boundaries, show- 
ing very good agreement across the entire range of T for 
b = 0.98. The spin correlations in phase C decay ac- 
cording to Eq. with r\ varying continuously between 
rj = 1/4 at T Ct i and 77 = 1/9 at T Ci2 H3- At the quantum 
critical point one has the 3D XY result 77 ps 0.040 |27| . 
The ordered phase (O) is characterized by finite sublat- 
tice magnetizations (V3/2, -V3/2, 0) which follows from 
<jjk = cos(QR jfe + 2n(ujk - uo)/3) = cos(QR jfc + 7r/6) 
for flat strings (ujk = 0) with uo = —a/4 in class (+ — 0). 
This kind of order is consistent with simulations [6j . 

An exciting perspective is that the string analogy ap- 
plies also to frustrated Ising systems with general cou- 



plings and quenched disorder examined in recent exper- 
iments [2^. A glassy state of the strings, similar to a 
vortex glass [2|J , would be an interesting possibility. 
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